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I. INTRODUCTION 



Daniele ColosEI and Robert Oecklll 
Instituto de Matemdticas, Universidad Nacional Autonoma de Mexico, 
Campus Morelia, CP. 58190, Morelia, Michoacdn, Mexico 

We investigate the question of unitarity of evolution between hypersurfaces in quantum field theory 
in curved spacetime from the perspective of the general boundary formulation. Unitarity thus means 
unitarity of the quantum operator that maps the state space associated with one hypersurface to 
the state space associated with the other hypersurface. Working in Klein-Gordon theory, we find 
that such an evolution is generically unitary given a one-to-one correspondence between classical 
solutions in neighborhoods of the respective hypersurfaces. This covers the case of pairs of Cauchy 
hypersurfaces, but also certain cases where hypersurfaces are timelike. The tools we use are the 
Schrodinger representation and the Feynman path integral. 

\ PACS numbers: 04.62.-|-v, ll.lO.-z 
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^sj ' Time evolution is implemented by unitary operators in quantum theory in order for a consistent probability in- 
terpretation to be possible, where probabilities are conserved in time. In non-relativistic quantum mechanics this is 
pH generally not an issue. Also in flat-spacetime quantum field theory, the unitarity of time evolution between equal-time 
[ hypersurfaces presents no problems. However, in quantum field theory it should also be possible to consider more 
^ • general spacelike hypersurfaces and have unitary evolution between them. In curved spacetime such a more general 
1^ I setting becomes a necessity. 

. The first explicit description of evolution between general spacelike hypersurfaces in quantum field theory was 
given by Tomonaga [1] and Schwinger |2|. They proposed a functional differential equation, known as the Tomonaga- 
, Schwinger equation that was to describe the infinitesimal evolution of wave functions between spacelike hypersurfaces. 
^ ' It has since been questioned not only whether such a heuristic equation can actually be made sense of, but also whether 
, a unitary evolution between spacelike hypersurfaces can be implemented in principle Q. 

I The purpose of the present article is to answer the latter question in the affirmative. What is more, we aim to 
' show that unitary evolution actually occurs generically, at least in the free Klein-Gordon theory (possibly with source 
^—i , terms). Moreover, unitary "evolution" also occurs generically between certain classes of timelike hypersurfaces. In 
' order to put this into perspective, it is necessary to clarify the conceptual framework on which we are building. 

We shall adopt here the point of view furnished by the general boundary formulation 0, which has recently 
emerged as a viable description of quantum fields offering a new perspective on their dynamics |6l-[l0|. In particular, 
this means that a priori there is not any single Hilbert space of states for the quantum field theory. Instead, a Hilbert 
space of states is associated with each hypersurface in spacetime. In the present paper these will be leaves of a 
certain foliation. The states of these Hilbert spaces may be thought of as encoding information about physics near the 
associated hypersurface. Probability conservation and unitarity takes on a more general meaning than is usual ITl| . 
d Indeed, whenever the dynamics of the theory links states on two hypersurfaces in a one-to-one correspondence, the 
associated quantum operator should be unitary. In a quantization context this comes from a one-to-one correspondence 
between classical solutions near the hypersurfaces. Such a correspondence occurs of course in the standard case of 
spacelike Cauchy hypersurfaces. However, it may also occur for certain families of timelike hypersurfaces. Indeed, we 
show in the present article that given such a correspondence on the classical level (and certain additional assumptions, 
in particular about the metric) the corresponding quantum operator relating the two hypersurfaces is unitary. In such 
a setting a single Hilbert space can be recovered a posteriori, if desired, by identifying via these unitary operators all 
Hilbert spaces associated to hypersurfaces that are in correspondence with a given one. 

The technical tools we shall use in this paper are quite simple and of limited sophistication. The point of view we are 
taking is naturally realized in the Schrodinger representation, where the state space associated with a hypersurface is 
simply the space of wave functions on field configurations on the hypersurface. Combining this with the Feynman path 
integral provides a convenient and sufficiently explicit description of quantum dynamics between hypersurfaces. In 
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particular, this allows to avoid any infinitesimal treatment of hypersurface deformations as in the Tomonaga-Schwinger 
approach, which is potentially problematic. 

We start in Section with a treatment of the classical Klein-Gordon equation, finding the action in the relevant 
spacetime region in terms of boundary data. We proceed in Section IIIII to lay out the framework of the quantum 
theory, starting with the free theory and its unitarity in Section [III Al Next follows a discussion of the compatibility 
of the vacuum wave function with the previous results in Section IIIIBI We then turn to a treatment of the free 
theory with a source term in Section [ill C[ showing how the results generalize to this case. We proceed with a short 
discussion of the theory with perturbative interactions in Section rill D[ without however trying to resolve the question 
of unitarity in this case. Finally in Section IIVI we illustrate our treatment with standard examples in Minkowski 
space. Throughout the paper, we assume for simplicity that certain fundamental operators commute, as happens 
in all examples that have been considered so far. The explicit form of some expressions in the more general case is 
contained in Appendix [XI In Appendix [B] we show how some of the fundamental operators are expressed in terms of 
Green's functions. 



II. CLASSICAL THEORY 



Consider the free theory of a Klein-Gordon field (f) with mass m propagating on a four-dimensional Lorentzian 
spacetime with line element given by 

ds^ = g^^dx^dx'', = 0,1,2,3. (1) 

The action in a spacetime region M is 



^ -Im 



where the integration is extended over the region M and we use the notation 9^ ~ d/dx^, and g = det g^,y. The 
equation of motion is obtained by varying the action with respect to the field and setting the variation equal to zero, 
yielding the Klein-Gordon equation 



[V\9\g''''d.)+mA<j,ix)=0. (3) 



The action for a solution of the Klein-Gordon equation can be computed performing an integration by parts in Q 
and using equation ([3]), yielding 



S 



where s denotes some coordinates on the boundary dM of M, g*^^^ is the determinant of the induced metric on the 
boundary and is the outward normal to dM . 

From now on we suppose a smooth coordinate system which defines a foliation covering all or part of space- 
time. Thus, X S M'^ are coordinates on the leaves while i is a variable indexing the leaves. We suppose that t E I, 
where / might be a finite interval, a semi-infinite interval or all of K. We identify x^ = t and x* = x^. We also require 
the metric to be block diagonal with respect to the foliation, i.e., g"* = = for all i G {1, 2, 3}. Note however, that 
in spite of the suggestive notation we do not necessarily require that be a temporal and x* be spatial coordinates. 

We consider a spacetime region M bounded by the disjoint union of two hypersurfaces that arise as leaves of the 
foHation, namely T, = {{t,x) ■ t = and T, = {{t,x) : t = £,}. We denote this spacetime region by The action 

evaluated on this region then takes the form, 

^lU] ^^^^IJ {\/\9^'^9nHi 3L) {dt<f) il x) - ^jW^m, X) {dtcty) (C, x)) , (5) 

where g and g denote the metric restricted to the hypersurfaces S and S respectively. We shall be interested in the 
value of S^^ |j as a function of boundary field configurations 



(6) 
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In order to obtain an explicit expression we will need to make further assumptions. In particular, we assume that 
solutions in J X K'^ give rise to a one-to-one correspondence between solutions in small neighborhoods of any two 
hypersurfaces with fixed t e /. In the case where the coordinate system gives rise to a foliation by spacelike hyper- 
surfaces this is the usual Cauchy property. Fixing a particular way to specify "initial" data r]a{x) and r]f)(x) we can 
write a solution of the Klein- Gordon as 

Ht,x) = {Xa{t)r,a){x) + iXbit)r,b){x), (7) 

where each Xi(t) is a linear operator from the space of initial data rji to solutions evaluated on hypersurfaces at fixed 
values of t. In the following we shall assume for simplicity that all these operators commute among each other and 
that certain operators satisfy a symmetry property. In all examples we have considered so far, this turns out to be the 
case. If the operators do not commute or the symmetry property is not satisfied, the computations become somewhat 
more involved, see Appendix [X] Inverting the matrix operator equation 



we obtain,^ 



(8) 



(9) 



where A(^, ^) Xa{(,)Xb{£^) — Xa{£,)Xi,{£^). Reinserting this into ([7]) yields the solution in / x as a function of the 
boundary data, 



This allows to evaluate ([S]), resulting in, 

S 



(</>) = ^ / d'sL 0) Q , (11) 



where the W^^ |j is a 2x2 matrix with elements W^^'^^^ , = 1, 2), given by 

= - J|g(3)g00| ^iM) ^(1,2) ^ _ /: (3)„00| Milll 



[€,€] V 1^ ^ I A(^,^) ^ A(^,^) 

where 

Ai(C,|) :=atA(t,e)|^^^ A2(e,e) :-9tA(e,0L=^- (13) 
The space of smooth solutions of ([3]) on M'' x / is equipped with the following symplectic form 

^^(0i>2) = ^ / d^x^\gi^)g^°\{cp,dtct>2-cl>2dt4>i). (14) 

which is independent of the choice of leaf S of the foliation. (See [l^l. In case of the leaves being spacelike, this is 
just the standard symplectic form.) This implies that the operator 

W y^|.g(3)g00| A2(t,<) = Ai(i,i) (15) 

is independent of t. Note that this immediately implies W^^ |j — W^^ . 



The boundary value problem we are considering here does of course in general not have a unique solution. Indeed, it would usually 
have infinitely many solutions. Thus, there are arbitrary choices involved in the inversion process we are describing. However, any such 
choice yields the same result in the moment we evaluate the action, i.e., once we arrive at ((TTJ. We therefore allow ourselves to gloss 
over this detail as we did on previous occasions 0, |^ when we showed that standard results are correctly reproduced. 
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III. QUANTUM THEORY AND UNITARITY 



The passage to the quantum theory is implemented by the Feynman path integral prescription. Moreover, the 
quantum dynamics of the field is described in the Schrodinger representation, where the quantum states are wave 
functionals on the space of field configurations. Thus, with a given spacetime hypersurface S we associate the space 
of state ■Hs of wave functions of field configurations on E. This state space carries the following inner product. 



/ I?^^eMV's(¥'), (16) 



where the integral is over all field configurations on the E. Amplitudes pu '■ T~LdM — >■ C are associated to spacetime 
regions M. State spaces and amplitudes satisfy a number of consistency conditions, see Q or Q. 

Consider now as above regions [^,^] that are bounded by hypersurfaces S at ^ and E at ^. In this case, the field 
propagator associated with the spacetime region is formally defined as 

%|](^,^)=/' l?0e'^i«^«i(^\ (17) 

where S^^ (0) is the action of the field in the region ^] and the integration is extended over all field configurations 

(j) that reduce to the boundary configurations ^p and (p on the boundary hypersurfaces E and E respectively. All 
the information on the dynamical evolution of the field between boundary configurations (p and is encoded in the 
propagator P7|) . The amplitude associated with the region and a state ® & H*^ = 'H^j^ is then 



PM{^i(^il^^) ^ J d(pd0ip({ip)ilj^{ip)Z^^^^^{ip,ip). (18) 

The main consistency condition for amplitudes is the composition property. In the present context this property 
means that the composition of two amplitudes for evolving in the t variable first from to ^2 and then from ^2 to ^3 
equals the amplitude for the direct evolution from ^1 to ^3. This translates into the following identity for propagators: 



%1,«3](</'1>'P3) = y X>(p2%i,C2](¥'l,^2)%2,«3](V'2,</'3)- (19) 

We recall from 0] or Q that unitarity of the quantum evolution between hypersurfaces E (at t = £_) and E (at 
t — ^) is equivalent to the equation 



dip I] {(fi, (p) (<^', (p) = 5{lp - tf'), (20) 

where Lp, tp' are field configurations on the hypersurface E and ip are field configurations on the hypersurface E that 
are integrated over. The delta function is meant with respect to the integral over field configurations on E. As is easy 
to see, this equation precisely guarantees that the inner product (jl6p of two states a,i t — remains the same when 
each states is evolved via the propagator Z^^ |j from t = ^ to t — 



A. Free Theory 

The first case we are interested in is the free theory with the action given by ^ . We shall use a subscript in the 
relevant quantities to indicate that we are referring to this case. The propagator ([TT]) can be evaluated by shifting the 
integration variable by a classical solution, 0ci, matching the boundary configurations in ^] = EU E, i.e. 0ci|s = <y5 
and 4>c\\± — 'P- Explicitly, 

where the normalization factor is formally given by 

iV[^_|]_o=/ X>0e'^[s-si-o^"^\ (22) 
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Explicit calculation shows that the composition property (|19p of the propagators (PT|) is satisfied if the following 
identity for the normalization factors (|22p holds, 



Using the identity 



we obtain 



Vip exp ( -i y d^a;(^A(^j = det 



^K-«3l.o%^,,,],o%^,C3l,o = det 



[?1,«3 



V 



27r 



iU/(1^2) ^(1,2) 

27r 27r 



This suggests the following solution for the normalization factor, which we shall henceforth adopt. 



(23) 



(24) 



(25) 



(26) 



We turn to the unitarity condition (PO)) . Since the action is real we have, 
Thus, 



N, 



lU],o 



(1,2), 



K,a,o 
5(93' - (p). 



det 



US 



2tt 



(27) 



where we used the normalization factor N^^ |j ^ as found in 



B. Compatibility of the vacuum 



In the following we are going to verify that a notion of vacuum, if it exists, is compatible with the structures we 
have found for the free theory, including the unitarity condition. According to the axioms of the general boundary 
formulation Q , a vacuum consists of a particular state associated with each leaf of the foliation. These vacuum states 
have to be related by evolution between leaves. In the free theory, this amounts to 

(28) 



(29) 



Making a Gaussian ansatz for the form of the vacuum state. 



V'c,o(¥') = Q exp 



1 



d?x'p{x){A^^){x) , 
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where is some operator and is a normalization factor ensuring that the vacuum state is normahzed, via 



ICfP =det 



2tt 



(30) 



Evolving the vacuum state from t — to t — ^ with the field propagator yields the consistency condition, 



^l = ^[«,l].o^« det 



A^ - iW 



(14) 



2tt 



-1/2 



(31) 



As we shall see, this equation can indeed be satisfied and allows us to fix relative phases for the normalization factors 
Ct ■ The general form of the operator At has been derived in [l^ , 



At^-iJ\g(^)g°° 



dticaXa{t)+cMt)) 

CaXa{t) + CbXb{t) ' 



(32) 



where Ca and Cf, are complex numbers, characterizing the vacuum wave function. In order for the vacuum state to be 
normalizable, these complex numbers must satisfy the condition c^Cb — cEca 0. One easily checks that 



Ct = det 



v/|g(^)g°°|i(c^Cb-c^cJAi(t,t) 

27r{CaXa{t) + CbXbit))^ 

provides a solution to both (|5D|) and pTjl . where equation (1^^ has been used. 



(33) 



C. Interaction with a source field 



We consider in this section the interaction with a real source field /i. The action takes the form 

SM,ti.{(t>) = S'a/,o(^) + / d^'x ^/\g\^i{x)4>{x). 



(34) 



M 



As before we are interested in the field propagator in the region [^,^] where we can write the action in terms of 
boundary field configurations as 



= "^Kil^oW + / d^^(A'?(£)¥'(£) + M|(£)^(^)) 



with 



LLclx) :— / dt V 9 



We can then write the corresponding propagator as 

N, 



The normalization factor can be expressed as 



d a; vlsl t^i^) oi{x) 



where a is a solution of the inhomogeneous Klein-Gordon equation 

(□ + TO^)a = II, 



(35) 

(36) 
(37) 

(38) 

(39) 
(40) 
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with boundary conditions, 



a can be written as 



WA(C,e) 



(41) 



(42) 



Given these ingredients one may exphcitly verify that psp satisfies the composition property. 
Checking the unitarity condition (|20p yields, 



d<^^[j_j-]_o(V5,<^)-^[j_|]_o(9'',^)exp ( i / d^^^5(^) (<y9'(£) - (^(2)) j , 



N, 



= <5(¥,-^') 



exp ( i / d^xfi^ix) {tp'{x) - <y9(x)) I S{ip - ip') 



(43) 
(44) 



In the last step we have used that the quotient (|39p has modulus one since the integrand on the right hand side of 
is real. This confirms unitarity in the presence of a source field. 



D. Perturbative interaction 



Finally, we consider the perturbatively interacting theory described by the action 



Sm,v{cI>) = Sm.o{4>) + / d^x.M^\V{x,ct>{x)), 

J M 

where V is an arbitrary potential. We may write this using functional derivatives as, 

exp(i5Af,v((/>)) = exp ^i j d^x^/\g(x)\ V ^x, ~i ^^'|_^^ ^ ^ exp (iS'M,^(0)) 

The propagator for the spacetime region ^] can then be written as 

5 



ti=0 



.[,,^.j_^(^,<^)=exp(i/^^^^_^dxy(^.,-i^ 



=^[5^|],o('^> 0) exp I i d xV \x, -i 



5 



N, 



(45) 



(46) 



A ^ ^ ; ; /\r exp 1 1 / d'^x^\g{x)\ n{x) (l){x) 



M=0 



^K,|],o('^' exp I i / . d 2^ ^ ( a;, -i 



X exp I 1 / d^x\p\g{x)\ix(x) (\)(x) 



5^{x) 



fj.=0 



exp 



US 



d'^x\/\g{x)\ fj,{x) Ufj, (x) 



(47) 



where in the last step the quotient of normalization factors p9p has been substituted, and the function is given by 
(j42p . where the subscript is meant to emphasize that depends on /i. The use of this functional derivative technique 
to express the field propagator in the presence of a perturbative interaction turns out to be inadequate to study the 
unitarity with the method proposed in this work, namely via the composition of the propagators (|20p . An alternative 
way to analyze the evolution for a general interacting theory would be to study the properties of the S-matrix. We 
shall elaborate on this elsewhere. 
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IV. EXAMPLES 



Although the present article is written very much in order to provide a context for unitary evolution in curved 
spacetime, we limit ourselves in this example section to fiat Minkowski spacetime. However, we emphasize instead 
the fact that the "evolution" does not have to be temporal. Thus, after considering the standard setting of evolution 
between equal-time hypersurfaces we switch to a unitary "evolution" in a radial direction as considered in ■ We 
do not explicitly repeat the calculations involved in showing the composition property, unitarity etc., but merely 
indicate the involved operators and normalization factors. 



A. Time-interval region in Minkowski space 

The foliation of Minkowski spacetime is the standard one here, indexed by a global time coordinate i S M, while 
the leaves of the foliation have coordinates ^ e K'^. The initial data hypersurface is located at i = 0. We exemplify 
two different choices for initial data. If we choose (77a, ?7h) to be the positive and negative energy components of the 
field at time t = we have 



Xa{t) = e-"^\ Xb{t) = e'"* with u ^-A^^ + m2, (48) 

Aj; being the Laplacian in the coordinates x. If we choose {rjaTTjij) instead to represent the initial value and its 
temporal derivative at time t = we have instead, 

X,(i)=cos(ic^i), Xb{t) = ^-^^^. (49) 
For a time interval [ii,i2] the matrix W\t^ t^^ of (fT^ results to be 

^ (cosu:[t,-H) -1 \ 

sinw(i2-ii)V ~^ cosu(t2 - h) ) ' ^ 

while the normalization factor N\^ti.t2] of (EH) turns out to be 



%i.t2]=det : — . (51) 

^ ^ \27r sma;(t2 - 



B. Hypercylinder in Minkowski space 

We consider now a different foliation of Minkowski space, see 0, Q for details. We use coordinates {t, r, f2) where 
t is the usual time coordinate, while r and 17 are spherical coordinates in space. (f2 is a collective coordinate for the 
angles on the 2-sphere.) The leaves of the foliation are now the timelike hypersurfaces of constant r, where r S (0, 00). 
(We are missing the time axis in this foliation.) Note that the Minkowski metric diagonalizes in the required way. In 
the present context it makes sense to divide the classical solutions into those that converge when r — >■ and those 
that diverge when r — >■ 0. Using rja to encode initial data of solutions that converge and ry^ for solutions that diverge, 
we can define 

Xa(r)e'^*yr(^^) ai{E,r)e'^'Yr{n), Xb(r) e'^'F^ (f^) hi{E,r)e'^'Yr{n), (52) 
where Y™ are the usual spherical harmonics, G M and a; and hi are defined as 



ji{r^/ E"^ - TV?), if E'^>w?, , , _ ( ni{rVE'^ - m'^), if E'^ > 



ai{E,r) = < ■-■1^2^ 2 and bi(E,r)^\._\ /-^ •-■1^2^ 2 (^3) 

where j;, nj, and are the spherical Bessel functions of the first and second kind and the modified spherical Bessel 
functions of the first and second kind respectively. Given < R < R, the matrix W^^ of (fT^ is the operator valued 
2x2 matrix 

i (-R^ai{E,R,R) 1/p 

Si{E,R,R) V 1/P ~R^ai{E,R,R) 



2 



^[RA-JJ^r-^A 1/J ' _n.2^.^^m)' (54) 
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where 



r VE^ - m2, if > m2, 
\Wm^-E^, if < m^, 



(55) 



and the functions 5i and ct; are to be understood as operators defined as 



Si {E, R, R) = ai {E, R)hi {E, R) - bi {E, R)ai {E, R), ai {E, R, R) = ai {E, R)b[ {E, R) - hi {E, R)a[ [E, R), (56) 



where the prime indicates derivative with respect to the second argument. The normahzation factor A^r^ ^, of (|26p is 



In the present paper we have investigated the quantum evolution of Schrodinger wave functions of Klein-Gordon 
theory along special foliations of (a part of) curved spacetime. We have worked out the field propagators between 
leaves of the foliation and found that they yield unitary "evolution" operators between the corresponding spaces of 
wave functions. The key required properties of the foliation were two: (a) There has to be a one-to-one correspondence 
between classical solutions in a neighborhood of any pair of leaves, (b) The spacetime metric has to take a block 
diagonal form with respect to the foliation. These conditions can be naturally realized for Cauchy hypersurfaces 
in globally hyperbolic manifolds. However, as we have shown, there are also interesting examples involving unitary 
"evolution" between timelike hypersurfaces. As already hinted at in the introduction, this implements probability 
conservation "in space" rather than in time, see Q. 

Although we have used foliations for convenience, it is clear from our construction that the amplitude (and thus 
the corresponding "evolution" operator) for any pair of hypersurfaces is manifestly independent of the rest of the 
foliation. The construction of the field propagator that determines the amplitude did not involve information about 
the intermediate foliation. Indeed, no such intermediate foliation needs to exist. (However, condition (a) above 
needs to be satisfied for the pair of hypersurfaces and condition (b) needs to be satisfied in neighborhoods of the 
hypersurfaces.) 

It is now time to clarify differences to other approaches at unitary evolution in curved spacetime that have run 
into difhculties. As already mentioned in the introduction the key difference to the original approach of Tomonaga 
and Schwinger is more of a technical than of a conceptual nature. The problems inherent in describing hypersurface 
deformations in terms of infinitesimal generators and the ambiguities associated to "integrating" these generators are 
simply avoided in a Feynman path integral approach. Instead of seeking a generalization of the Schrodinger equation 
we directly look for a generalization of transition amplitudes. 

Another approach that has been followed in the literature is to try to generalize the fact that the Poincare group 
acts on the state space of quantum field theory in Minkowski space. Of course, a generic Lorentzian manifold will 
not have any isometrics. However, supposing there is a Killing vector field leaving the Lagrangian invariant, this 
will induce an action on the phase space of the classical field theory by symplectic transformations. One may then 
hope to implement such an action on a Hilbert space for the quantum field theory and interpret it as the generator 
of evolutions along the vector field [TBL Indeed, one may even try actions coming from suitable more general 
(non-isometric) spacetime transformations 3]. In either case, serious difficulties or even no-go theorems have resulted. 
Let us emphasize that in our approach no such action is obtained. Indeed, the only natural way to relate the Hilbert 
spaces associated with different hypersurfaces in our approach is precisely via the unitary operators representing 
evolution between them. Absent any other way to relate the Hilbert spaces, no non-trivial action can be constructed. 

A deficiency of our paper may be seen in the fact that the concrete examples we have presented concern only 
Minkowski spacetime, rather than the curved spacetimes our treatment is mainly aimed at. However, some curved 
spacetime examples have already been partially worked out in the present setting (isl . [l6j and others are to follow. 



In the main body of the article, we have assumed that all operators of the form Xa{t), Xii{t') commute and that 

(i i) 

the operators satisfy a symmetry condition of the form 




(57) 



V. 



CONCLUSIONS 



Appendix A: Non-commuting operators 




(AI) 
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where (p, tp' are field configuration data. If this is not the case, the computations become somewhat more involved. In 
the foUowing we wiU present merely the resulting elements of the matrix W^^ |j that reduce in the commutative case 

to the expressions given in (fT^ . The general result can be written as follows: 

3| (^X',X-'D-^XtX-' - X',X-'D-^XaX-^) , (A2) 

~\{X',X-'D-^-X',X-^D-^), (A3) 

(X',X-^D-^X,X^^ - XiX^'D-'XaX-') , (A4) 

(X^X^'D-' KX-'D-') . (A5) 

Here we have used the abbreviated notation Xi Xi{(^) and Xi :— Xi{S^) and D :— X^X^^ — XaX^^. What we still 
assume is that Xa{t) and Xi,{t) are invertible. However, this is merely a consequence of the one-to-one correspondence 
between initial data and solutions near any hypersurface of the foliation. We also assume that D is invertible without 
being able to offer here an a priori justification. 



Appendix B: Using Green's Functions 



Various operators appearing in our treatment can be made more explicit by using specific Green's functions of the 
Klein-Gordon equation. 

The inverse of the symplectic structure is a function G : (/ x K.-^) x (/ x M'') M which is antisymmetric in 
its two arguments and satisfies the condition, 

2ni<f>,Gy)^^{y), where Gyix) -.^ G{x,y), (Bl) 

and where (j) is an arbitrary classical solution (see [HI)- Indeed, this implies that G is actually a Green's function 
and solves the Klein-Gordon equation in both of its arguments. Furthermore, it is easy to see that it satisfies the 
boundary conditions 

G{{t,x),{t,y)) = V2,yeR^V^e/, (B2) 
^ 15^3) {t, 2L)9°° {t, x) I a,G((r, x), {t, y))\r=t = -^^'^ (s - y) yx,yem',tel. (B3) 

Consider now the operators Xa and Xi, of Section [ill As a first step we express them via integral kernels Xa and 
Xb, i.e., 

(X,(i)r;0(£) - J d\X,{t,x,iMy). (B4) 

In order to relate these to a Green's function we have to specify what the "initial" data r/a, rjb actually represent. We 
consider the following choice here: rja specifies the value of the field on the hypersurface t = to while ry^ specifies the 
normal derivative of the field on the same hypersurface: 

Va{x) ^ (f>{to,x), (B5) 



Vbix) = V\9"°{to,x)\drcbiT,x) . (B6) 

T = to 

We obtain 

Xa{t,x,y) = -J\g(^Hto,y)g°°{to,y)\drG{{T,y),{t,xj) , (B7) 

- V - - - r=to 



Xb(t,x,y) = ^|.g(3)(to,y)| G{{to,y), {t,x)), (B8) 
where G is precisely the Green's function considered above. 
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The operators reconstructing a classical solution from its values on two hypersurfaces t = ^ and t = £ of the foliation 
as in equation (|1U|) can also be related directly to a certain Green's function, different from G. (See [ij], where such 
methods were used in a related context.) As a first step we write these operators in terms of integral kernels, 



( mt) 



(p] ix)= I d^yK^^^^^{t,x,y)(p{y), 



(B9) 
(BlO) 



Now consider the symmetric Green's function Gj^ |j : (/ x M^) x (/ x M"^) — > M with the following boundary conditions: 



G[^,|] ((£,£), = V3i,y e R^V^ e /, 

The integral kernels K and L can then be written as 



(Bll) 
(B12) 

(B13) 
(B14) 



The operators can also be written as integral kernels in a nicely symmetric way. Define 

{w[l'2^) [x) = j <^yW^^{x,y)m. 
Then, combining equations p2p with the expressions above yields, 



W, 



(B15) 



(B16) 



Here we have set £i := £ and £2 := The relative sign when i ^ j originates from the opposite orientation of the two 
hypersurfaces as boundaries of the intermediate region. 
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